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Abstract 

Accurate  information  on  the  temperature  field  and  associated  heat  transfer  rates  are  particularly  important  in  devising  appropriate  heat  and  water 
management  strategies  in  proton  exchange  membrane  (PEM)  fuel  cells.  An  important  parameter  in  fuel  cell  performance  analysis  is  the  effective 
thermal  conductivity  of  the  gas  diffusion  layer  (GDL).  Estimation  of  the  effective  thermal  conductivity  is  complicated  because  of  the  random  nature 
of  the  GDL  micro  structure.  In  the  present  study,  a  compact  analytical  model  for  evaluating  the  effective  thermal  conductivity  of  fibrous  GDLs 
is  developed.  The  model  accounts  for  conduction  in  both  the  solid  fibrous  matrix  and  in  the  gas  phase;  the  spreading  resistance  associated  with 
the  contact  area  between  overlapping  fibers;  gas  rarefaction  effects  in  microgaps;  and  salient  geometric  and  mechanical  features  including  fiber 
orientation  and  compressive  forces  due  to  cell/stack  clamping.  The  model  predictions  are  in  good  agreement  with  existing  experimental  data  over 
a  wide  range  of  porosities.  Parametric  studies  are  performed  using  the  proposed  model  to  investigate  the  effect  of  bipolar  plate  pressure,  aspect 
ratio,  fiber  diameter,  fiber  angle,  and  operating  temperature. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Electrochemical  energy  conversion  in  hydrogen  fuel  cells  is 
an  exothermic  process  that  results  in  significant  temperature 
variations  [1,2],  Accurate  information  on  the  temperature  field 
and  associated  heat  transfer  rates  are  particularly  important  in 
devising  appropriate  heat  and  water  management  strategies  in 
proton  exchange  membrane  (PEM)  fuel  cells,  as  the  temper¬ 
ature  field  affects  relative  humidity,  membrane  water  content, 
and  reaction  kinetics,  as  well  as  durability.  One  of  the  main 
fuel  cell  components  in  this  respect  is  the  porous  gas  trans¬ 
port  layer,  commonly  referred  to  the  gas  diffusion  layer  (GDL). 
GDLs  employed  in  PEM  fuel  cells  typically  consist  of  a  fibrous 
structure  in  the  form  of  a  thin  “paper”  or  “woven  cloth”,  see 
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Fig.  1.  The  GDL  provides  five  key  functions  for  a  PEM  fuel 
cell;  (1)  mechanical  support,  (2)  electronic  conductivity,  (3)  heat 
removal,  (4)  reactant  access  to  catalyst  layers,  and  (5)  product 
removal  [3]. 

The  porous  nature  of  GDL  micro  structure,  makes  it  necessary 
to  define  an  effective  thermal  conductivity,  a  transport  parame¬ 
ter  that  plays  an  important  role  in  fuel  cell  performance  analysis 
[4]  and  that  is  required  in  computational  models  [5].  In  addi¬ 
tion  to  being  porous,  GDLs  are  anisotropic,  which  adds  to  the 
complexity  of  characterizing  the  effective  thermal  conductivity. 

Ramousse  et  al.  [4]  recently  investigated  the  effective  thermal 
conductivity  of  non-woven  carbon  felt  GDLs  and  estimated  the 
conductivity  bounds  using  a  model  connecting  the  two  phases 
(solid  and  gas)  in  series  or  parallel.  The  model  as  well  as  their 
experimental  measurements  yielded  conductivity  values  that  are 
at  least  one  order  of  magnitude  lower  than  most  values  reported 
in  the  literature.  Ramousse  et  al.  [4]  also  noted  that  due  to  contact 
and  constriction  resistances  between  carbon  fibers,  the  effective 
thermal  conductivity  of  carbon  felts  are  much  lower  than  pure 
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Nomenclature 

A  area  (m2) 

a ,  b  major  and  minor  semi  axes  of  elliptical  contact 

area  (m) 

d  fiber  diameter  (m) 

E  Young’s  modulus  (Pa) 

E'  effective  elastic  modulus  (Pa) 

F  contact  load  (N) 

F\  integral  function  of  (p'  p"~l),  Eq.  (8) 

GDL  gas  diffusion  layer 

K(-)  complete  elliptic  integral  of  the  first  kind 

k  thermal  conductivity  ( Wm- 1  K“ 1 ) 

keff  effective  thermal  conductivity  (Wm" 1  K“ 1 ) 
keff0  effective  thermal  conductivity  of  the  reference 
basic  cell  (Wm-1  K-1) 

k*  non-dimensional  effective  thermal  conductivity, 

fctff&effo 

Z  distance  between  two  adjacent  fibers  in  re¬ 

direction  (Fig.  3)  (m) 

/jbp  bipolar  pressure  (Pa) 

Pg  gas  pressure  (Pa) 

Pqdl  GDL  pressure  (Pa) 

Pr  Prandtl  number  (-) 

Qgc  heat  transfer  rate  through  gas  filled  gap  (W) 

R  thermal  resistance  (KW_1) 

Rco  constriction  resistance  (KW_1) 

Rsp  spreading  resistance  (KW_1) 

T  temperature  (K) 

Vs  fiber  (solid)  volume  of  basic  cell  (m3) 

Vtot  total  volume  of  basic  cell  (m3) 

w  distance  between  two  adjacent  fibers  in  the  y- 
direction  (Fig.  3)  (m) 

Greek 

a  thermal  accommodation  parameter 

P  fluid  property  parameter,  Eq.  (17) 

S(x)  local  gap  thickness  (m) 

e  porosity  (-) 

rj  modulus  of  elliptic  integral  (-) 

Y  heat  capacity  ratio  (-) 

A  mean  free  path  of  gas  molecules  (m) 

A.  ratio  of  relative  radii  of  curvature  (p'  p"-1) 

p,  ratio  of  molecular  weights  of  the  gas  and  the  solid, 

Mg  M~l 

9  angle  between  two  fibers  (rad) 

p' ,  p"  major  and  minor  relative  radii  of  curvature  (m) 

p[ ,  p2  principal  radii  of  curvature  (m) 

pe  equivalent  radius  of  curvature  of  the  contacting 
surfaces  (m) 

v  Poisson’s  ratio  (-) 

§  aspect  ratio  (wl“ 1 ) 

Subscripts 

0  reference  state 

1  bottom  block  of  the  basic  cell 


2  top  block  of  the  basic  cell 

oo  standard  condition  state 

c  contact  plane 

g  gas 

gc  gas  filled  gap 

max  maximum  value 

s  solid  (carbon  fiber) 

t  upper  boundary  of  the  top  block 

tot  total  value 


carbon,  and  used  Danes  and  Bardon  [6]  correlation  to  estimate 
the  effective  thermal  conductivity  of  the  solid  phase. 

Khandelwal  and  Mench  [7]  measured  the  through-plane  ther¬ 
mal  conductivity  of  GDLs.  They  examined  two  different  types 
of  commercial  GDLs  with  a  variety  of  thickness  and  porosity. 
They  studied  the  effect  of  temperature  and  polytetrafluoro  ethy¬ 
lene  (PTFE)  content  on  the  effective  thermal  conductivity,  and 
obtained  values  in  close  agreement  with  the  manufacturer  data 
sheet.  The  experimental  data  reported  in  the  literature  for  effec¬ 
tive  thermal  conductivity  spans  a  wide  range  of  values,  0.1-1. 6 
Wm-1  K_1,  and  there  is  clearly  need  to  better  understand  of 
the  possible  sources  of  inconsistency  [1,4], 

Our  literature  review  indicates  the  need  for  a  general  model 
that  can  accurately  predict  the  effective  thermal  conductivity  of 
GDL,  and  its  trends  as  parameters  varied,  since  no  reliable  corre¬ 
lations  are  available  and  there  is  lack  of  data  and  understanding 
on  the  effect  of  geometric  parameters  such  as  tortuosity,  radius 
of  contact  area  between  fibers,  and  angle  between  fibers.  The 
objectives  of  the  present  work  are: 

•  Develop  and  verify  a  comprehensive  analytical  model  that  can 
predict  the  effective  thermal  conductivity  of  GDLs  and  that 
captures  accurately  the  trends  observed  in  experimental  data. 

•  Investigate  the  effect  of  relevant  geometrical,  thermal,  and 
mechanical  parameters  involved  and  identify  the  controlling 
parameters. 


Fig.  1.  SEM  image  of  a  Toray  GDL. 
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Following  the  approach  used  successfully  in  several  applica¬ 
tions  such  as  spherical  packed  beds  by  Bahrami  et  al.  [8] ,  a  “basic 
cell”  is  taken  to  represent  the  fibrous  media,  i.e.  the  structure  is 
assumed  to  be  repeated  throughout  the  GDL.  Each  cell  is  made 
up  of  contact  regions.  A  contact  region  is  composed  of  a  contact 
area  between  two  portions  of  fibers,  surrounded  by  a  gas  (air) 
layer.  A  thermal  resistance  model  is  then  constructed  taking  into 
account  the  basic  conduction  processes  through  both  the  solid 
fibrous  matrix  and  the  gas  phase  as  well  as  important  phenom¬ 
ena  including  spreading  resistance  associated  with  the  contact 
area  between  overlapping  fibers  and  gas  rarefaction  effects  in 
microgaps. 

The  basic  cell  approach  breaks  the  problem  into  distinct  con¬ 
duction  paths,  the  contact  between  two  fibers,  the  gas  layer 
between  fibers;  and  calculates  the  conductivity  of  the  medium 
as  a  series/parallel  combination  of  the  individual  resistances  for 
those  paths.  The  advantage  of  this  approach  is  that  it  readily 
allows  evaluation  of  the  relative  contributions  of  each  conduc¬ 
tion  path  as  a  function  of  the  medium  properties  [8]. 

The  scheme  of  the  present  approach  to  evaluate  the  effec¬ 
tive  thermal  conductivity  is  shown  in  Fig.  2.  The  first  step  in 
estimating  the  effective  thermal  conductivity  is  the  reconstruc¬ 
tion  of  the  GDL  geometrical  structure.  The  GDL  is  represented 
as  cylindrical  carbon  fibers  that  are  equally  spaced  horizon¬ 
tally  and  stacked  vertically  to  form  mechanical  contacts.  Fig.  3. 
The  next  step  is  mechanical  modeling  of  the  contacting  fibers. 
The  Hertzian  theory  [9]  is  used  to  evaluate  the  contact  area 
between  fibers,  and  a  thermal  resistance  network  is  constructed 
to  account  for  the  effective  thermal  conductivity,  allowing  ana¬ 
lytic  determination  of  the  effective  thermal  conductivity.  The 
results  of  the  model  are  compared  to  experimental  data.  More¬ 
over,  parametric  study  is  then  performed  to  investigate  the 
effect  of  key  parameters  on  the  effective  thermal  conductivity  of 
GDLs. 

2.  Model  development 

Both  electrical  and  thermal  conductivity  of  carbon  paper 
GDLs  are  orthotropic  [4,10],  with  in-plane  conductivity  that  are 


Fig.  2.  Model  development. 


an  order  of  magnitude  higher  than  the  through  plane  value.  The 
thermal  field  and  heat  transfer  rates  depend  on  a  variety  of  factors 
including,  geometry,  material  properties  of  the  various  compo¬ 
nents  and  operating  conditions,  the  heat  transfer  in  the  GDL 
is  however  generally  limited  by  the  through  plane  conductiv¬ 
ity  value  on  which  we  focus  our  analysis.  The  model  considers 
the  GDL  to  consist  of  a  periodic  fibrous  micro  structure  and 
assumes: 


(1)  3-D  repeating  basic  cell,  Fig.  3. 

(2)  Steady  state  one-dimensional  heat  transfer. 

(3)  Negligible  natural  convection;  justified  by  the  Grashof  num¬ 
ber  for  a  typical  GDL  with  fiber  diameter  of  8.5  pm  which 
is  in  order  of  10-6  and  is  significantly  lower  than  2500,  the 
limit  for  natural  convection  [11], 

(4)  No  radiation  heat  transfer.  PEM  fuel  cells  typically  operate 
between  60  and  90  °C,  and  the  contribution  of  radiation  is 
small  and  can  be  neglected. 


Fig.  3.  (a)  Front  view  and  (b)  top  view  of  the  geometrical  model  of  GDL. 
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(5)  Regular  fiber  surfaces  (no  roughness  or  out-of-flatness  for 
contacting  fibers). 

Based  on  these  assumptions,  we  propose  a  resistance  network 
model  for  conduction  through  the  solid  and  gas  phases  which 
accounts  for  geometric  structure,  effect  of  compression,  gas  rar¬ 
efaction  effect  (in  microgaps  between  fibers),  and  spreading 
resistance. 

2.1.  Geometrical  model 


Fig.  4.  Pressure  distribution  on  the  bipolar  plate  and  GDL. 


Fig.  1  shows  SEM  image  of  a  Toray  GDL  that  clearly 
illustrates  the  random  and  anisotropic  structure.  The  proposed 
geometrical  model  is  an  idealization  shown  in  Fig.  3  and  consists 
of  uniformly  sized  equally  spaced  cylindrical  fibers  immersed 
in  stagnant  air.  The  fibers  angle,  0,  is  variable  in  this  model. 

The  porosity  is  defined  as: 


(1) 


where  Vs  and  Vtot  are  the  volume  of  fibers  and  the  basic  cell, 
respectively.  Calculating  these  volumes  based  on  the  basic  cell 
geometry  in  Fig.  3  yields: 


jtd  [l  +  (10/ cos  9)' 
T  [  Iw  . 


(2) 


2.2.  Mechanical  model 


Thermal  energy  transfers  from  one  fiber  to  another  through 
the  contact  interface,  and  resistance  to  heat  conduction  depends 
on  the  contact  area  dimensions.  In  order  to  determine  the  contact 
area  dimensions,  the  Hertzian  theory  [9]  is  used  in  the  present 
study. 

The  general  shape  of  the  contact  area  is  elliptical;  when  9  =  0, 
the  contact  area  becomes  circular.  Applying  the  Hertzian  theory 
[9],  the  semi-axes  of  the  contact  area  are  given  by: 


where  a  and  b  are  the  major  and  minor  semi-axes  of  the  ellip¬ 
tical  contact  area,  respectively;  pe  is  the  equivalent  radius  of 
curvature,  pe  =  ~J  P'  P"  [9];  and  E'  is  a  modulus  incorporating 
the  fibers  Young’s  modulus  and  Poisson’s  ratio. 


(5) 


p'  and  p"  are  the  major  and  minor  relative  radii  of  curvature  at 
the  contact  point  expressed  as  [9]: 


„  _  d 

P  ~  V2(l  —  cos  20)  +  2 


(6) 


P  (4/d)  -  (1/p")  (?) 

F\ ,  the  parameter  used  in  Eq.  (3),  is  a  complex  integral  function 
of  (X  =  p'/p ")  [9].  We  correlate  this  integral  and  propose  the 
following  relationship: 


Fi 


19.1-v/I 

1  +  16.76-s/k  +  1.34A. 


(8) 


The  accuracy  of  the  relationship  is  within  0.08%. 

In  order  to  determine  the  contact  area  dimensions,  the  mag¬ 
nitude  of  the  contact  force  is  required.  This  force  F  can  be 
evaluated  from  the  clamping  pressure  applied  via  the  fuel  cell 
bipolar  plates.  The  land/flow  channel  area  ratio  used  in  PEM 
fuel  cells  is  optimized  to  balance  electrical  conduction  and  mass 
transport  and  is  typically  of  order  1  as  shown  in  Fig.  4.  Thus  the 
maximum  pressure  to  which  the  GDL  is  subjected  is  twice  that 
of  the  bipolar  plate,  i.e.  Pgdl  =  2/jbp.  As  shown  in  Fig.  3,  a 
cell  with  the  cross-sectional  area  of  Iw  consists  of  four  contact 
points;  therefore,  the  corresponding  maximum  force  on  each 
contact  is: 


Fm  ax 


/JGDL  lw 

4 


(9) 


2.3.  Thermal  model 


Based  on  the  assumptions  discussed  in  Section  2,  a  thermal 
resistance  network  corresponding  to  the  basic  cell  is  constructed 
as  shown  in  Fig.  5  by  considering  the  top  and  bottom  blocks  of 
the  basic  cell  structure.  The  thermal  resistances  in  the  network 
represent  the  heat  transfer  paths  through  the  gas  and  carbon 
fibers,  see  Fig.  6.  The  solid  bulk  resistance  is  small  compared 
with  the  gap  resistance  and  its  effect  on  the  total  resistance  is 
negligible  and  not  accounted. 

2.3.1.  Thermal  constriction/spreading  resistance 

Thermal  constriction/spreading  resistance  is  defined  as  the 
difference  between  the  average  temperature  of  the  contact  area 
and  the  average  temperature  of  the  heat  source/sink,  which  is 
located  far  from  the  contact  area,  divided  by  the  total  heat  flow 
rate  Q  [12], 

AT 

Rco  =  Rsp  =  -  (10) 

If  the  contact  areas  are  small  compared  with  the  distance 
separating  them,  the  heat  source  on  a  half-space  solution  can 


yyyyyyyyyyyy. 


204 


E.  Sadeghi  et  al.  /Journal  of  Power  Sources  179  (2008)  200-208 


Fig.  5.  Thermal  resistance  network  for  the  top  and  bottom  blocks  of  the  basic 
cell. 

be  used  [13].  Fig.  7  illustrates  the  geometry  of  a  circular  heat 
source  on  a  half-space. 

The  spreading  resistance  for  an  isothermal  elliptical  contact 
area  can  be  determined  analytically  [14]: 

**=7^*™  (U) 


where  rj  =  \J  1  —  ( b/a )2  and  K{rj)  is  the  complete  elliptic  inte¬ 
gral  of  the  first  kind  of  modulus  17. 

f71!2  d  t 

K(jf)  =  /  ,  „  ,  :  (12) 

Jo  v  [1  —  t]2  siir  t] 

2.3.2.  Gas  resistance 

Gas  resistance  can  be  decomposed  into  two  parallel  resis¬ 
tances  for  each  block,  see  Fig.  6.  Kennard  [15]  modeled  the  gas 
conduction  between  two  parallel  plates  for  temperature  jump 


as  <2gc  =  kgAAT/(8  +  M).  Yovanovich  [16]  showed  that  this 
expression  can  be  used  for  all  possible  regimes.  We  assume  that 
heat  conduction  at  each  differential  element,  dx,  is  similar  to  that 
between  parallel  plates;  therefore,  using  Kennard’s  expression 
for  each  differential  element,  we  find: 


d(2gc  —  kg  (  2  ) 


A  T{x) 
S(x)  +  M 


dx 


(13) 


The  gas  parameter  M  is  defined 


M  =  up  A 


(14) 


The  mean  free  path  A  of  the  gas  molecules  can  be  expressed 
in  terms  of  /lgj00,  the  mean  free  path  of  gas  molecules  at  the 
reference  state. 


with  Tg  oo  =  o  =  l  atm.  The  thermal  accom¬ 
modation  par  ie  fluid  property  parameter  P  are 

defined  by: 


(17) 


where  y  is  the  ratio  of  specific  heats,  Pr  is  the  Prandtl  number, 
and  a\,  U2  are  thermal  accommodation  coefficients  of  the  top 
and  bottom  surfaces.  Flere,  the  top  surface  corresponds  to  carbon 
fiber  and  the  bottom  to  the  gas,  i.e.  u\  =  as  and  «2  =  1 , 


/  2  -  as ) 


+  1 


(18) 


Song  and  Yovanovich  [17]  experimentally  correlated  the  ther¬ 
mal  accommodation  coefficient. 


«s  =  exp  —0.57 


/rs-273\l/  1.4Mg  \ 

V  273  )\  \  6.8  +  1  AMg  ) 


/Ts-273 
V  273 


)]} 


2.4/r, 

(1  +  A*-)2 


Fig.  6.  Thermal  resistances  for  the  top  block  of  the  basic  cell. 
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where  ji  =  Mg/Ms;  Mg  and  Ms  are  molecular  weights  of  the 
gas  and  the  solid. 

The  total  heat  flow  through  the  gas-filled  gap,  gas  filled 
between  the  quarter  cylindrical  fiber  and  the  separating  plane,  is 
given  by  the  integral 


AT(x) 
S(x)  +  M 


(20) 


The  thermal  resistance  of  the  gas-filled  gap,  Rgc,2,  is  defined 
in  terms  of  the  temperature  difference  between  two  bounding 
surfaces,  A  Tc. 


Qgc 

1  A  Tc  '' 


The  local  gap  thickness  can  be  defined  based  on  the  geometry 
of  the  contact  interface. 


<$(*)  «  \/x  ~(x  cos  Of 


(22) 


Table  1 

Properties  of  air 


^g,oo  Gun) 

£g,oo  (W/mK) 

7g  (K) 

Pg  (kPa)  Pr  y 

0.07  [20] 

0.03  [4] 

350  [21] 

101.3  [4]  0.7164  1.398 

Table  2 

Carbon  fiber  properties 

d  (|xm) 

ks  (W/mK) 

v 

E  (GPa)  Pbp  (kPa) 

8.5 

120  [4] 

0.3  [22] 

210  [23]  482 

program  are  given  in  Tables  1  and  2.  The  gas  phase  is  taken 
as  air  to  correspond  to  available  experiments  used  for  valida¬ 
tion. 

3.1.  Model  validation 


Considering  isothermal  bounding  surfaces  ,  AT {x)  =  ATC, 
the  thermal  resistance,  Eq.  (21)  reduces  to: 


1  rd/2cose  dx 

Rgc,2  ~kg\  2J  JQ  8{x)  +  M 


(23) 


The  thermal  resistance  Rgc, 2  can  then  be  calculated  by  sub¬ 
stituting  Eqs.  (14)  and  (22)  into  Eq.  (23).  This  equation  can  also 
be  used  for  Rgp_  by  setting  8{x)  =  d/2  and  a  =  2  (both  surfaces 
are  gas): 


J_  =  ,  (v>\  ((/-(<// cos  0))/2) 
% 2  g  V  2  )  {d/2)  +  a/)A 


(24) 


The  thermal  resistances  for  the  bottom  block  can  be  obtained 
following  the  same  procedure. 


/ 1  cos  0  \  fd/2  cos0  dx 

V  2  )  Jo  8{x)+M 

(25) 

' l  cos 0\  ((to- </)/2  cos 0) 

,  2  )  {d/2)  +  a/3 A 

(26) 

2.3.3.  Effective  thermal  conductivity 

Once  the  individual  resistances  are  determined,  the  thermal 
resistance  network  shown  in  Fig. 5  is  used  to  evaluate  the  total 
resistance  of  the  basic  cell: 


[  1  ,  1  1  I 

r'j 

r 1  +  1  + 1 1 

k  p  Rgc, 2  Rg,i\ 

1  +l 

[fico  figc,l  fig,lj 

(27) 


The  effective  thermal  conductivity  of  the  GDL  is  given  by: 


,  d  4  d 

eff  ~~  fitotA  ~~  lwRm 


(28) 


3.  Results  and  discussion 


The  model  was  implemented  into  a  Mathematica  script  [18] 
to  allow  convenient  parametric  studies  and  analysis.  The  ther¬ 
mophysical  properties  of  the  gas  and  carbon  fibers  used  in  the 


Fig.  8  compares  the  model  to  experimental  data  from  a  variety 
of  sources  obtained  over  a  range  of  porosities.  As  shown  in  Fig.  8, 
there  is  good  agreement  between  the  model  and  experimental 
data  with  an  average  difference  of  approximately  7.5%  when 
0  =  0.  The  model  results  for  three  arbitrarily  chosen  angles  are 
also  shown  in  Fig.  8.  The  case  0  =  0  yields  better  overall  agree¬ 
ment  with  experimental  data  and  for  small  fiber  angles,  the  effect 
of  angle  is  negligible.  In  an  actual  GDL,  the  fiber  angle  distri¬ 
bution  is  random  and  it  is  expected  that  the  average  value  would 
correspond  to  the  orthogonal  arrangement,  0  =  0,  as  the  present 
model  suggests.  This  is  corroborated  by  the  recent  results  of 
Van  Doormal  [19]  who  performed  Lattice  Boltzmann  simula¬ 
tion  in  reconstructed  GDLs  and  observed  that  the  orthogonal 
arrangement  yields  permeability  values  that  are  indistinguish¬ 
able  from  those  computed  using  a  random  fiber  arrangement. 
The  case  0  =  0  is  therefore  selected  as  the  reference  case  for  the 
parametric  studies. 


Fig.  8.  Comparison  of  predicted  and  measured  effective  conductivities  over  a 
range  of  GDL  porosities.  (See  Refs.  [4,7,24,25]). 
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e (deg) 


Fig.  9.  Effect  of  fibers  angle  ( 0 )  on  effective  thermal  conductivity. 

3.2.  Parametric  study 

The  proposed  model  can  be  conveniently  used  to  systemati¬ 
cally  investigate  the  influence  of  key  GDL  parameters/properties 
on  its  effective  thermal  conductivity.  Parametric  studies  were 
conducted  by  varying  fibers  angle,  bipolar  pressure,  aspect  ratio, 
fiber  diameter,  and  operation  temperature.  When  a  parameter  is 
studied,  the  rest  of  mentioned  parameters  are  kept  constant. 

3.2.1.  Fibers  angle 

The  preferred  conduction  path  corresponds  to  smallest  resis¬ 
tance  which  is  through  the  contact  between  fibers.  This  path 
includes  spreading  resistance;  and  is  hence  expected  to  have 
a  significant  influence  on  the  effective  thermal  conductivity. 
The  effect  of  fibers  angle  on  the  non-dimensional  properties 
of  the  basic  cell  is  shown  in  Fig.  9.  The  properties  are  non- 
dimensionalized  with  respect  to  the  reference  case,  see  Table  3. 
As  shown  in  Fig.  9,  the  contact  area  increases  with  increasing 
0  in  turn  results  in  a  reduction  in  spreading  and  consequently, 
total  resistances.  However  this  is  counterbalanced  by  an  increase 
in  the  basic  cell  area,  and  because  the  latter  dominates  over  the 
resistance  in  this  case,  the  net  effect  is  a  reduction  of  the  effective 
thermal  conductivity.  This  effect  becomes  more  pronounced  for 
higher  fiber  angles. 

3.2.2.  Bipolar  pressure 

The  effect  of  bipolar  pressure  on  the  thermal  conductivity  of 
GDL  is  shown  in  Fig.  10.  Higher  bipolar  plate  pressures  result 
in  higher  thermal  conductivities.  A  higher  pressure  leads  to  an 

Table  3 

The  reference  case  parameters 

0  w  =  l  (pm)  dofiun)  £  Pbp  (kPa)  Tg  (K) 

0  33.38  8.5  0.8  482  350 


Fig.  10.  Effect  of  bipolar  plate  pressure  on  effective  thermal  conductivity. 

increase  in  the  contact  area  that  in  turn  leads  to  a  decrease  in  the 
spreading  resistance,  and  thus  higher  effective  thermal  conduc¬ 
tivity.  It  should  be  noted  that  as  a  result  of  increasing  pressure, 
the  height  of  the  basic  cell  is  expected  to  decreases  due  to  elas¬ 
tic  compression.  However,  this  effect  is  small,  i.e.  a  reduction  of 
less  than  1%  in  fiber  diameter,  thus  is  neglected  in  this  study.  In 
an  operating  fuel  cell,  the  compressive  force  to  which  the  GDL 
is  subjected  is  expected  to  vary  from  a  maximum  under  the  land 
area  of  the  bipolar  plate  to  a  minimum  under  the  centre  of  the 
flow  channel.  The  sensitivity  to  pressure  shown  by  the  model 
suggests  that  the  effective  conductivity  is  non-homogeneous  and 
this  should  be  accounted  for  in  comprehensive  fuel  cell  models. 

3.2.3.  Aspect  ratio 

Fig.  1 1  shows  the  effect  of  aspect  ratio  £  =  w/ 1  for  a  GDL 
with  a  porosity  e  =  0.8.  As  shown,  the  lower  the  aspect  ratio, 
the  lower  the  effective  thermal  conductivity.  When  the  aspect 
ratio  is  reduced,  in  order  to  maintain  the  same  porosity,  Z  has  to 
be  increased  while  w  is  decreased,  see  Eq.  (2).  This  leads  to  a 
larger  basic  cell  area,  and  since  the  bipolar  pressure  is  kept  con¬ 
stant  results  in  larger  contact  force  and  hence  a  lower  spreading 
resistance.  Thus  the  total  resistance  decreases  when  the  porosity 
is  constant,  but  this  is  counteracted  by  a  proportionally  larger 
increase  in  the  basic  cell  area,  and  thus  a  lower  effective  thermal 
conductivity. 

3.2.4.  Fiber  diameter 

The  effect  of  carbon  fiber  diameter  on  the  effective  thermal 
conductivity  at  a  constant  porosity  is  shown  in  Fig.  12.  The  total 
resistance  decreases  with  an  increase  in  the  fiber  diameter.  To 
keep  the  porosity  constant,  however,  the  length  and  width  of  the 
basic  cell  have  to  be  increased  as  the  diameter  increases.  Con¬ 
sequently,  the  area  of  the  basic  cell  becomes  larger  and  the  total 
thermal  resistance  decreases.  However,  a  larger  diameter  leads 
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Fig.  1 1 .  Effect  of  aspect  ratio  on  effective  thermal  conductivity. 

to  a  larger  basic  cell  area  which  negatively  impacts  the  effec¬ 
tive  thermal  conductivity.  Because  the  effect  of  the  cell  height 
extension  and  total  resistance  reduction  is  higher  than  the  area 
enlargement,  the  net  effect  is  a  bit  increase  in  the  effective  ther¬ 
mal  conductivity  with  fiber  diameter.  In  the  typical  range  of 
5-10  p.m  for  the  fiber  diameter,  the  effective  thermal  conductiv¬ 
ity  remains  approximately  constant. 

3.2.5.  Operating  temperature 

Fig.  13  shows  the  effect  of  temperature  on  the  effective 
thermal  conductivity.  The  typical  operating  temperatures  of 
automotive  PEM  fuel  cell  is  in  the  80-90  °  C  range  while  air 
cooled  and  air  breathing  cells  operate  at  lower  temperatures; 


d'  =  d/d0 

Fig.  12.  Effect  of  fiber  diameter  on  effective  thermal  conductivity. 


T  (C) 

Fig.  13.  Effect  of  the  operating  temperature  on  effective  thermal  conductivity. 


we  thus  investigate  variations  in  the  15-100°  C  range.  The 
thermal  and  mechanical  properties  of  the  fiber  are  assumed  to 
remain  constant  within  this  range;  however,  the  effects  of  tem¬ 
perature  variations  are  considered  in  gas(air)  thermophysical 
properties  e.g.  M,  k,  etc.  As  shown  in  Fig.  13,  the  gap  resistance 
RgC  increases  slightly  with  temperature,  whereas  the  gas  resis¬ 
tance  Rg  decreases.  The  two  effects  balance  each  other  and  also 
the  spreading  resistance  does  not  vary,  therefore,  the  effective 
thermal  conductivity  remains  approximately  constant. 

4.  Summary  and  conclusions 

A  compact  analytical  model  for  evaluating  the  effective  ther¬ 
mal  conductivity  of  fibrous  GDLs  has  been  developed.  The 
model  accounts  for  the  salient  geometric  features,  the  effect 
of  mechanical  compression,  and  spreading  resistance  through 
fibers.  The  model  predictions  are  in  good  agreement  with  exist¬ 
ing  experimental  data  over  a  wide  range  of  porosities,  0.3  < 
e  <  0.8.  Parametric  studies  have  been  performed  using  the  pro¬ 
posed  model  to  investigate  the  trends  and  effects  of  bipolar  plate 
pressure,  aspect  ratio,  fiber  diameter,  fiber  angle,  and  operating 
temperature.  The  highlights  of  the  analysis  are: 

•  Constriction/spreading  resistance  is  the  controlling  compo¬ 
nent  of  the  total  resistance. 

•  Orthogonal  arrangement  of  fibers,  0  =  0  yields  better  overall 
agreement  with  experimental  data  and  corroborated  by  the 
recent  results  of  Van  Doormal  [19]. 

•  The  influence  of  fiber  angle  0  on  the  effective  thermal  con¬ 
ductivity  decreases  at  higher  porosities. 

•  Higher  bipolar  pressure  significantly  improves  the  effective 
thermal  conductivity. 

•  Reducing  the  aspect  ratio,  §  =  w/l,  to  approximately  0.7  has  a 
negligible  impact  on  the  effective  thermal  conductivity.  How¬ 
ever,  for  f  <  0.3,  the  effect  of  aspect  ratio  becomes  important. 
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The  analysis  indicates  that  the  best  effective  thermal  conduc¬ 
tivity  is  achieved  when  £  =  1  (square  basic  cell). 

•  Neither  changes  in  fiber  diameter  (5-10pun)  nor  operating 
temperature  for  the  range  of  15-100  °C  have  any  significant 
effect  on  the  effective  thermal  conductivity. 

The  compact  model  presented  here,  reproduces  faithfully  the 
effects  of  many  operational  and  geometrical  parameters  on  effec¬ 
tive  thermal  conductivity.  The  model  can  be  used  to  guide  the 
design  of  improved  GDLs,  and  can  be  readily  implemented  into 
fuel  cell  models  that  require  specification  of  the  effective  thermal 
conductivity. 

Acknowledgements 

The  authors  are  grateful  for  the  financial  support  of  the  Nat¬ 
ural  Sciences  and  Engineering  Research  Council  (NSERC)  of 
Canada,  and  the  Canada  Research  Chairs  Program.  The  SEM 
image  in  Fig.  1  was  provided  by  Aimy  Bazylak. 

References 

[1]  N.  Djilali,  D.  Lu,  Int.  J.  Therm.  Sci.  41  (2002)  29-40. 

[2]  A.  Hakenjos,  H.  Muenter,  U.  Wittstadt,  C.  Hebling,  J.  Power  Sources  131 
(2004)  213-216. 

[3]  J.P.  Feser,  A.K.  Prasad,  S.G.  Advani,  J.  Power  Sources  162  (2006) 
1226-1231. 

[4]  J.  Ramousse.S.  Didierjean,  P.  Lottin,  D.  Maillet,  Int.  J.  Therm.  Sci.  47 
(2008)  1-8. 

[5]  B.R.  Sivertsen,  N.  Djilali,  J.  Power  Sources  141  (1)  (2005)  65-78. 


[6]  F.  Danes,  J.P.  Bardon,  Revue  Generate  de  Thermique  36  (1997)  302-311. 

[7]  M.  Khandelwal,  M.M.  Mench,  J.  Power  Sources  161  (2006)  1106-1115. 

[8]  M.  Bahrami,  M.M.  Yovanovich,  J.R.  Culham,  Int.  J.  Heat  Mass  Transfer 
49 (2006)  3691-3701. 

[9]  K.L.  Johnson,  Contact  Mechanics,  Cambridge  University  Press,  London, 
UK,  1985  (Chapter  4). 

[10]  W.  Sun,  B.A.  Peppley,  K.  Karan,  J.  Power  Sources  144  (1)  (2005)  42-53. 

[11]  V.S.  Arpaci,  P.S.  Larsen,  Convection  Heat  Transfer,  Prentice-Hall,  Engle¬ 
wood  Cliffs,  NJ,  1984  (Chapter  4). 

[12]  H.S.  Carslaw,  J.C.  Jaeger,  Conduction  of  Heat  in  Solids,  second  ed.,  Oxford 
University  Press,  London,  UK,  1959. 

[13]  A.M.  Clausing,  B.T.  Chao,  J.  Heat  Transfer  87  (1965)  243-251. 

[14]  M.M.  Yovanovich,  E.E.  Marotta,  A.  Bejian,  D.  Kraus,  in:  Heat  Transfer 
Hand  Book,  John  Wiley  and  Sons  Inc,  New  York,  2003  (Chapter  4). 

[15]  E.H.  Kennard,  Kinetic  Theory  of  Gases,  McGraw-Hill,  New  York,  1938 
(Chapter  8). 

[16]  M.M.  Yovanovich,  Spacecraft  Radiative  Transfer  and  Temperature  Control, 
AIAA,  New  York,  83  (1992)  83-95. 

[17]  S.  Song,  M.M.  Yovanovich,  Fundamentals  of  Conduction  and  Recent 
Developments  in  Contact  Resistance  ASME  HTD,  69  (1987)  107-1 15. 

[18]  Wolfram  Mathematica  6,  Wolfram  Research,  Inc.,  1988-2007. 

[19]  M.  Van  Doormaal,  Determination  of  permeability  in  fibrous  porous  media 
using  the  lattice  Boltzmann  method  with  application  to  PEM  fuel  cells,  M. 
Eng.  thesis,  Queen’s  University,  Kingston,  2006. 

[20]  M.V.  Williams,  E.  Begg,  L.  Bonville,  H.R.  Kunz,  J.M.  Fenton,  J.  Elec- 
trochem.  Soc.,  A  151  (8)  (2004)  1173-1180. 

[21]  K.  Rajashekara,  J.  IEEE  41  (3)  (2005). 

[22]  H.  Miyagawa,  C.  Sato,  T.  Masea,  E.  Drowna,  J.  Mater.  Sci.  Eng.,  A  412 
(2005)  88-92. 

[23]  A.  Kusoglu,  A.M.  Karlsson,  M.H.  Santare,  S.  Cleghom,  W.B.  Johnson,  J. 
Power  Sources  161  (2006)  987-996. 

[24]  A.  Rowe,  X.  Li,  J.  Power  Sources  102  (2001)  82-96. 

[25]  H.  Ju,  C.Y.  Wang,  S.  Cleghom,  U.  Beuscherb,  J.  Electrochem.  Soc.,  A  52 
(2005) 1645-1653. 


